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Abstract 

We have successfully grown sizable single crystals of La2Cui_ z (Zn,Mg) z 04 with 
up to nearly half of the magnetic Cu sites replaced by non-magnetic Zn and Mg. Neu- 
tron scattering, SQUID magnetometry, and complementary quantum Monte Carlo 
(QMC) simulations demonstrate that this material is an excellent model system for 
the study of site percolation of the square-lattice Heisenberg antiferromagnet (SL- 
HAF) in the quantum-spin limit S = 1/2. Carefully oxygen-reduced samples exhibit 
Neel order up to the percolation threshold for site dilution, z p « 40.7%. For z > 10%, 
the material exhibits a low-temperature tetragonal (LTT) structural phase, with a 
transition temperature that increases linearly with doping. Above z ~ 25%, Neel 
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order occurs in the LTT phase. Up to at least z = 35%, the Neel temperature 
Tn(z) of the experimental system corresponds to the temperature at which QMC 
indicates that the spin correlations for the nearest-neighbor S = 1/2 SLHAF have 
grown to approximately 100 lattice constants. Neutron scattering measurements of 
the static structure factor in the paramagnetic regime allow the determination of the 
two-dimensional spin correlations, which are found to be in excellent quantitative 
agreement with QMC over a wide common temperature and doping range. Neutron 
scattering and QMC results for the temperature dependence of the static structure 
factor amplitude S(tt,tt) are in good agreement as well. As the concentration of 
non-magnetic sites is increased, the magnetic correlation length £(T, z) crosses over 
from an exponential dependence on p s /T to power-law behavior in the temperature 
regime studied. Fits to a heuristic crossover form for £(T, z) allow an estimate of the 
spin stiffness, p s = p s (z), which approaches zero at z = z p . The combined experi- 
mental and numerical data presented here provide valuable quantitative information 
for tests of theories of the randomly-diluted 5=1/2 SLHAF. 

Key words: Magnetic dilution, two-dimensional, Heisenberg antiferromagnet, 
quantum phase transition, neutron scattering, quantum Monte Carlo, percolation 
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The discovery of high-temperature superconductivity in charge-carrier doped 
La2Cu04 has stimulated an enormous interest in the properties of low-dimensional 
quantum magnets. The main structural building block of La 2 Cu04 and of re- 
lated Mott insulators is a square lattice of Cu 2+ spin-1/2 magnetic ions that 
experience a strong nearest-neighbor (NN) antiferromagnetic superexchange 
(J 1550 K) mediated by intervening oxygens. Apart from small correction 
terms (such as anisotropies and a three-dimensional coupling), which lead to 
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Neel order at a non-zero temperature, the magnetic degrees of freedom of these 
Cu0 2 sheets are well-described by the square-lattice Heisenberg Hamiltonian 

H=j£S,-S„ (1) 

(id) 

where the sum is over NN sites, J is the antiferromagnetic Cu-O-Cu superex- 
change, and Sj is the S — 1/2 operator at the site i. 

Extensive experimental [1,2,3,4,5,6,7,8,9,10], quantum Monte Carlo [11,12,13], 
and theoretical [14,15,16,17,18,19,20,21,22,23,24] efforts have led to a good 
understanding of the quantum many-body physics described by Eq. 1. The 
discrete Hamiltonian Eq. (1) has been mapped to a continuum quantum non- 
linear a model (QNLcrM), and it has been established that the NN square- 
lattice Heisenberg antiferromagnet (SLHAF) exhibits a broken-symmetry ground 
state (and renormalized classical behavior) rather than a quantum-disordered 
ground state even in the extreme quantum-spin limit of spin-1/2 [14]. Numeri- 
cal and theoretical studies indicate that the quantum critical point separating 
these two ground states may be traversed, for example, upon the introduction 
of frustrating next-NN interactions [14,25] or by coupling two SLHAFs to form 
a bilayer [26]. Although the full spin Hamiltonian of La2Cu04 is believed to 
contain a frustrating next-NN exchange of about 0.05-0. 10 J [7,8,27], this is 
well below the value of 0.24 J needed to disorder the system [25]. 

Quantum phase transitions in the presence of disorder are the subject of con- 
siderable current interest. Magnetic systems are often of particular value in 
this context, as the combined effects of quantum fluctuations and quenched 
disorder can be studied with various analytical and numerical methods. Fur- 
thermore, model magnets, if they can be found in nature, can serve as an 
important testing ground of theoretical predictions. A significant amount of 
theoretical and numerical effort has been devoted to the relatively simple one- 
dimensional [28,29] and two-dimensional [30,31] Ising model in a transverse 
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field. For the 0(3) symmetric S — 1/2 SLHAF, the effects of a single impu- 
rity are well understood [32,33,34], but there are few theoretical results for 
finite impurity concentrations [35]. There have been suggestions that random 
site [30,36,37] and bond [38] dilution of the NN S = 1/2 SLHAF may lead 
to a non-trivial quantum phase transition. Extensive experimental results ex- 
ist for S = 5/2 up to very high impurity concentrations [39], but results for 
the extreme quantum-spin limit of S=l/2 have been limited to lower con- 
centrations [5,40,41,42,43,44,45,46,47,48]. In studies of site-diluted S = 5/2 
Heisenberg and Ising antiferromagnets [39,49], long-range order was found to 
disappear only above the percolation threshold z p ks 40.725% [50,51]. Earlier 
results for the Neel temperature of randomly-diluted La 2 Cu0 4 extrapolate 
to zero temperature at concentrations well below the percolation threshold 
[5,40,41,42,43,44,45,46,47]. These results, along with theoretical predictions 
[30,36,37] and numerical studies [52], suggested the possible existence of a new 
quantum critical point at z s =i/2 < z p . However, recent Monte Carlo simula- 
tions [53,54,55,56,57,58,59] indicate that the site diluted NN S = 1/2 SLHAF 
remains ordered up to the percolation threshold. While the zero-temperature 
transition appears to be a classical percolation transition, the intermediate- 
temperature properties of the spin- 1/2 system at z = z p seem to be controlled 
by the effective proximity to a new multicritical point [58,59] 

Very recently, we have been able to show that La 2 Cui_ z (Zn,Mg) z 4 is a model 
material for the study of the effects of nonmagnetic impurities in the S — 1/2 
SLHAF up to and beyond the site percolation threshold [56]. Figure 1 is 
a schematic of the randomly diluted Cu02 sheets in La2Cui_ z (Zn,Mg) z 04, 
showing the geometric transition from an infinite connected cluster to finite 
disconnected clusters as the percolation threshold is traversed. We have also 
performed quantum Monte Carlo (QMC) simulations of the randomly diluted 
NN S = 1/2 SLHAF, 
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n = j Y,PiPA-Sj, (2) 

where pi = 1 (pi = 0) on magnetic (nonmagnetic) sites. The experimen- 
tal results for the static structure factor in the 2D correlated parramag- 
netic region agree remarkably well with our numerics [56], indicating that 
La 2 Cui_ z (Zn,Mg) z 04 is indeed a very good model system. The present paper 
outlines some aspects of these results in more detail. 

0.1 Sample Growth and Characterization 

A lack of stoichiometric samples at high non-magnetic concentrations had 
limited earlier experimental work on randomly diluted L^CuOzi. In previous 
studies of polycrystalline samples the concentration of non-magnetic ions was 
as high as 25% [47], significantly below the percolation threshold. Single crystal 
results had been limited even further, to z ~ 15% and below [5,44,46]. Fur- 
thermore, the excess oxygen typically found in as-grown samples introduces 
holes into the copper-oxygen sheets which frustrate the antiferromagnetism 
and quickly destroy magnetic order [60]. Differing values for T N {z) indicate 
that this problem had not been fully resolved [5,40,41,42,43,44,45,46,47]. 

We have successfully grown single crystals of La2Cui_ z (Zn,Mg) z 04 by the opti- 
cal traveling-solvent floating-zone method. Zn 2+ and Mg 2+ are non-magnetic 
and have larger and smaller ionic radii, respectively, than Cu 2+ . By jointly 
substituting Zn and Mg for Cu, we succeeded in growing samples up to and 
beyond the percolation threshold [56]. The Zn content of our samples is ap- 
proximately 10%, while the Mg content varies. Typical single-grain sections 
are about 4 mm in diameter and 40 mm long. The composition of the samples 
was measured with electron probe microanalysis on sections cut from the ends 
of the crystals. Some crystals have small Cu/Zn/Mg concentration gradients 
along the length (Az = 1-2%), but all crystals are very uniform radially. The 
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mosaic widths are very good, 15' full width at half maximum (FWHM) or 
less, as measured by neutron diffraction. Samples were carefully annealed for 
24 hours at 900 to 950 °C in an Ar atmosphere to remove excess oxygen. Us- 
ing SQUID magnetometry, we checked that subsequent anneals do not further 
raise T/v, confirming that the first anneal was successful. 

The neutron scattering measurements presented here were performed at the 
National Institure of Standards and Technology (NIST) Center for Neutron 
Research (NCNR), using the BT2 and BT9 thermal instruments and the 
SPINS cold neutron instrument. 

0.2 Phase Diagram 

Near the Neel transition, the uniform susceptibility of La2Cu 1 _ z (Zn,Mg) z C , 4 be- 
comes coupled to the staggered susceptibility through the antisymmetric 
Dzyaloshinskii-Moriya term (not included in Eqs. 1 and 2), and the uniform 
susceptibility perpendicular to the Cu0 2 sheets increases to a cusp at T N 
[46,61]. In order to determine the Neel temperature, small pieces, a few mm 
on each side, were cut from the larger crystals for magnetometry. Figure 2a 
shows the uniform susceptibility for several samples measured with a SQUID 
magnetometer. Polycrystalline samples, with concentrations assumed to be 
equal to their nominal values, yielded similar results. At higher dilution lev- 
els, a large Curie-like component emerges [46], and the cusp at T N becomes 
more difficult to distinguish. 

Neutron diffraction measurements allow us to follow T^(z) at higher concen- 
trations. The intensity at a magnetic Bragg reflection is proportional to the 
square of the staggered magnetic moment, and the doping dependence of the 
extracted ordered moment [56] agrees rather well with recent theoretical [35] 
and numerical [57] results for the NN SLHAF. Figure 2b shows the tempera- 
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ture dependence of the magnetic scattering at the (1,0,0) magnetic reflection 
for several samples (we use orthorhombic (Bmab) notation throughout). The 
neutron diffraction results for T N agree with SQUID measurements at low 
and intermediate concentrations and indicate that Neel order persists up to 
at least z = 39%. 

We typically observe an additional temperature-independent non-magnetic 
signal which was subtracted from the data in Fig. 2b. This remnant signal 
is a combination of double scattering from the (1,0,2) and (0,0,2) nuclear 
Bragg reflections and a small A/2 contamination due to nuclear scattering of 
higher-energy neutrons. The double-scattering component is strongly energy- 
dependent, and decreases at lower neutron energies. Neutrons with half the 
wavelength (and four times the energy) of the primary beam will also be 
diffracted by the monochromator. These neutrons are filtered out with py- 
rolytic graphite on thermal instruments or berylium on SPINS, but small 
numbers may still pass through if the filter is not sufficiently thick, giving 
(2,0,0) structural Bragg scattering at the same angular position as the pri- 
mary beam at the (1,0,0) magnetic peak. 

The correction terms to Eqs. 1 and 2 are believed to include a small frus- 
trating next-NN exchange [7,8,27], which might in principle become more 
important at higher concentrations, where T/v is relatively small, and lead 
to low-temperature spin-glass behavior. The associated spin freezing would be 
expected to lead to a momentum broadening of the magnetic peaks. Moreover, 
the observed transition temperature would depend on the probe frequency. We 
have tested for the possibility of low-temperature spin-glass physics, and found 
in all samples that showed magnetic order (up to z — 39%), that the width 
of the observed magnetic peaks remained resolution limited. Furthermore, we 
measured the temperature dependence of the magnetic scattering in a 35% 
diluted sample using incident neutron energies of 3.5, 5, and 13 meV, with 
FWHM energy resolutions ranging from 0.08 to 0.72 meV, as shown in Fig. 3. 
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Since the transition temperature is independent of the neutron energy reso- 
lution over a relatively wide range of energies, we conclude that the observed 
behavior is a genuine Neel transition [62]. 

At high temperatures, pure La2Cu04 is in the so-called high-temperature 
tetragonal (HTT) phase (IA/mmm). Below about 530 K, the oxygen octa- 
hedra surrounding the copper ions tilt in a staggered fashion, creating an 
orthorhombic distortion of the CuC>2 planes and forming the low-temperature 
orthorhombic (LTO) phase (Bmab). In a 8% Zn-doped sample, we find that 
the HTT to LTO transition temperature increases to 577 K, but we have 
not followed this transition to higher concentrations, although previous ex- 
periments have found that it increases monotonically with both Zn and Mg 
doping [41]. Above z = 10%, we discovered a second structural transition 
into a low-temperature tetragonal (LTT) phase (PA 2 /ncm). Figure 4 shows 
longitudinal scans through the (2,0,0) position measured by neutron diffrac- 
tion above, at, and below the LTT transition temperature Tltt — Q3(5)K for 
z = 19%. The two peaks observed at higher temperature correspond to the 
(2,0,0) and (0,2,0) reflections. Both are observed in the same scan because 
of the presence of twin domains with different orthorhombic distortion direc- 
tions. In the LTT phase, only one reflection is observed. The redistribution of 
intensity from the orthorhombic peaks to the tetragonal peak without a shift 
in the peak positions indicates that this is a first-order transition. We note, 
that the LTT phase has previously been observed in (La,Ba) 2 Cu0 4 [63] and 
(La,Sr,Nd) 2 Cu0 4 [64]. 

The magnetic and structural phase diagram obtained from SQUID magne- 
tometry and neutron diffraction is shown in Fig. 5. Below z ~ 20%, our data 
agree well with several previous results [5,45,46,47]. Above this concentration, 
we find that Tn(z) deviates from a linear behavior, approaching zero only at 
the percolation threshold. Although random dilution weakens the tendency to 
order [56], it appears that quantum fluctuations are not strong enough to shift 
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the critical point: -25=1/2 = z p , within the uncertainty of our experiment. 

The pure two-dimensional Heisenberg antiferromagnet described by Eq. 1 can 
not exhibit long-range order at nonzero temperature. However, weak inter- 
plane couplings and anisotropies in La2Cui_ z (Zn,Mg) z 04 lead to Neel order 
at nonzero temperature approximately when a e ff £f D ~ 1, where £20 is the 
two-dimensional magnetic correlation length corresponding to Eq. 1 and a e ff 
is a suitable combination of the correction terms in the full spin Hamiltonian 
[7]. In pure L^CuOzi, this occurs when £ 2 .d ~ 100a, where a is the planar 
lattice constant [7]. Using QMC results for Eq. 2 (discussed in more detail 
below), we find remarkable agreement between the £,2d/ci = 100 contour and 
Tn(z) up to at least 35%. This is demonstrated in Fig. 5a. The continuous 
line indicates direct numerical results for Eq. 2, while the dashed line results 
from an extrapolation of numerical data at higher temperatures. The full spin 
Hamiltonian describing La2Cui_ z (Zn,Mg) z 04 should depend on the details of 
the crystal structure. Above z ~ 25%, Neel order occurs in the LTT phase. 
Nevertheless, T N {z) evolves smoothly with doping, and corresponds to £20/° = 
100 from QMC for Eq. 2 even when T N < T LT t, as demonstrated in Fig. 5. 
Consequently, any changes with doping in the full spin Hamiltonian must be 
very subtle. 

A nonzero next-NN exchange in La2Cui_ z (Zn,Mg) z 04 could shift the perco- 
lation threshold from the NN-only value of z p « 40.7%. Since the next-NN 
coupling in La 2 Cui_ z (Zn,Mg) z 4 is frustrating [8], it could, in principle, shift 
the critical point noticeably to a value below the percolation threshold (and 
lead to spin-glass physics, as discussed above) [39]. However, the value of the 
next-NN exchange is relatively small (0.05-0.10J) [8]. For the pure Hamilto- 
nian Eq. 1, an additional frustrating next-NN coupling of ps 0.24J is needed to 
disorder the ground state [25]. We note that recent spin- wave measurements at 
the magnetic zone boundary suggest that the dominant further-neighbor in- 
teraction is not a next-NN exchange but a four-spin ring exchange [10], which 



9 



does not extend connectivity beyond the NN percolation threshold. 
0.3 Static Structure Factor Measurement 

We have systematically studied the static structure factor in the paramagnetic 
phase of La2Cui„ z (Zn,Mg) z 04, which allowed us to determine the instanta- 
neous two-dimensional (2D) spin-spin correlation length £(z,T). Above T N , 
the 3D magnetic Bragg peaks become rods of 2D scattering. The equal-time 
structure factor S(q2D), where q2D is the 2D momentum transfer component in 
the CuC>2 sheets relative to the 2D magnetic zone center, was measured with 
neutron scattering in two-axis, energy-integrating scans across these rods [7]. 
Figure 6 shows representative data for two different samples. The measured 
peaks broaden as £ decreases, both with increasing temperature and increasing 
dilution. Correlation lengths were obtained from fits to a single 2D Lorentzian 

convoluted with the instrumental resolution. We note that well above T N) the 
spin system is effectively isotropic. However, near TV the anisotropy and 3D 
coupling terms in the full spin Hamiltonian lead to crossover physics. Our 
two-axis, energy-integrating experiment simultaneously measures both the in- 
plane and out-of-plane components of the static structure factor. Since our 
data do not allow for a line-shape analysis, we carried out fits to a single 2D 
Lorentzian. 

The correlation lengths are plotted versus J/T in Fig. 7a, and the static struc- 
ture factor amplitude S(n, it) is shown in Fig. 7b. Temperature is scaled by J 
= 135 meV, the antiferromagnetic superexchange energy of the pure system 
[3,4]. The data are cut off above T N by one standard deviation 4 K), as 
obtained from fits of the order parameter (Fig. 2b). From Fig. 7a it can be seen 
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that dilution significantly decreases the rate at which correlations grow as the 
system is cooled. At high concentrations, T) crosses over from exponential 
to power-law behavior. We note that the z = 40(2)% and 43(2)% samples do 
not exhibit Neel order at 1.4K (J/T ~ 1100), and the spin correlations appear 
to approach a constant zero-temperature value, as expected for z > z p . 

0.4 Quantum Monte Carlo 

In order to test the degree to which La2Cui_ z (Zn,Mg) z 04 is described by the 
randomly-diluted NN Heisenberg Hamiltonian we have performed quantum 
Monte Carlo (QMC) simulations to calculate £(z,T) for Eq. 2. We used a 
loop-cluster algorithm [56,65,66,67], with lattice sizes 10 to 20 times larger 
than the correlation length in order to avoid finite-size effects. We were able 
to simulate very large lattices of up to 1700x1700 sites, and to reach temper- 
atures as low as T = J/100. Previous simulations were confined to very small 
lattices and high temperatures [52,68]. At each temperature and concentra- 
tion, between 5 and 200 random configurations were averaged, with 10 4 to 10 5 
measurements per configuration. The QMC results for £(z,T), shown as black 
symbols in Fig. 7a, extend to higher temperatures and complement the exper- 
imental data. We find excellent quantitative agreement between the two up 
to the percolation threshold. We emphasize that this comparison contains no 
adjustable parameters, since J and a are known. Above the percolation thresh- 
old, the effective concentration of approximately 46% is slightly higher than 
the actual experimental value of 43(2)%. A possible origin for this might be 
a stronger relative influence of the next-NN term in the full spin Hamiltonian 
for z > z p . 

In Fig. 7b, the static structure factor amplitude S(n, it) from QMC (filled 
symbols) is shown together with the experimental results. Unlike the correla- 
tion length, the measured amplitude also depends on experimental conditions, 
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such as the effective illuminated sample volume and the neutron flux. Since 
the absolute value of S(tt, it) could not be determined with good accuracy, we 
have normalized the experimental data for each sample to match the respec- 
tive numerical values. Note that the temperature scale is not adjustable. The 
temperature dependence of S(ir, ir) can therefore still be compared between 
experiment and QMC, and we find good agreement. 



0.5 Theory 

The ground state of the pure NN SLHAF is ordered, but quantum fluctuations 
renormalize the spin- wave velocity, c = 2\^2SZ C (S) Ja, and spin-stiffness, p s = 
S 2 Z P (S)J, from their classical values (using units in which g/iB — ks — h — 1). 
For S = 1/2, the quantum renormalization factors Z c and Z p are known from 
theoretical and numerical studies, and £(z = 0,T) is given by [14,15] 



£ e c/a 



a 8 2ixp s 



i-i ^\ + o 



2 \2itp s ) \2nps 



(4) 



with c = 1. 657 J a and p s = 0.18 J [13]. Even though Eq. 4 is strictly valid 
only at asymptotically low temperatures [13,14,15,24], it agrees remarkably 
well with experiment [6,7,8] and numerics [12,13] in the range 2 < £/a < 200 
shown in Fig. 7a. 

The derivation of Eq. 4 involves mapping the discrete Hamiltonian (Eq. 1) 
to the continuum quantum nonlinear sigma model (QNLctM), and it assumes 
the existence of an ordered ground state and translational invariance. Ran- 
dom dilution breaks translational invariance and leads to defect rods in the 
Euclidian-time direction of the effective QNLaM. These rods break Lorentz 
invariance, and a QNLaM description may no longer be valid [14,35]. Further- 
more, we can see from Fig. 7a that at higher concentrations the temperature 
dependence of £(z,T) approaches power-law behavior. The modified form of 
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Eq. 4 



a 



e c/a e ^ps/T 
8 27rp s l + (47rp s /T) 



(5) 



with v T — 1 has been suggested for disorder-free systems approaching a quan- 
tum critical point [69]. This crossover formula interpolates between Eq. 4 at 
T < p s and £ ~ 1/T as p s approaches zero. Although recent studies indicate 
that z = z p is not a quantum critical point [57,58,59], random dilution reduces 
the spin stiffness and may be viewed as bringing the system closer to such a 
point in an extended parameter space. 

In a classical picture, at T = 0, z = z p is a (geometric and thermal) multicrit- 
ical point, and p s = as well as power-law behavior of £(z p ,T) are expected. 
For the S = 5/2 system Rb 2 (Mn,Mg)F 4 at and above the percolation thresh- 
old, the correlation length is well described by the form 



where £o(z) is the saturated zero-temperature length and £t(T) ~ T~ UT [70]. 
At the percolation threshold, the thermal exponent of this S = 5/2 system was 
found to be vr = 0.90(5). A model based on the growth of ID spin correlations 
along self-avoiding walks on the percolation cluster has been used successfully 
to describe these results for Rb 2 (Mn,Mg)F 4 [70] without adjustable parame- 
ters. However, this model assumes classical spin chains, and the fact that we 
obtain a different exponent for S — 1/2 (see below) is therefore not surprising. 

We have tested the extent to which Eqs. 4 and 5 describe our results by 
fitting the numerical data using c(z) and p s (z) as adjustable parameters. We 
find that fits to Eq. 4 give a good description of £(z,T), especially at low 
concentrations, and even at z — 31% for £/a > 8. However, the modified 
crossover form Eq. 5 even captures the high-temperature power-law behavior 
at higher concentrations, as shown in Fig. 7a. Results for c(z) and p s (z) from 
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fits of Eq. 5, including dilution levels not shown in Fig. 7, are shown in Fig. 
8. Using Eq. 4 results in large uncertainties for c(z), but both c(z) and p s {z) 
obtained using the two different forms agree within the errors. For the pure 
system, fits of QMC results [13] below £/a = 200 yield 27rp s (0) = 1.18(1) J 
and c(0) = 1.33(3) Ja, about 4% higher and 20% lower, respectively, than the 
most accurate estimate [13]. For z > z p , the combined QMC and experimental 
data were fit to Eq. 6. As can be seen in Fig. 7a, our results are described well 
by this form, and we obtain v T = 0.72(7) for z = 41%. 

An attempt to combine percolation theory with the QNLaM model predicts 
for the spin stiffness [37] 

9(0) 



Ps(0) 



A(z) 



1 



(7) 



where g(0) = 0.685 is the coupling constant corresponding to the NN S=l/2 
SLHAF at z = [14], A(z) is the bond dilution factor, and Poo(z) is the 
probability of a site belonging to the infinite cluster of spins. A(z) is well 
described up to z « 37% by A(z) « 1 - ixz + irz 2 /2 [71]. Up to z « 20%, 
-Poo (z) 7n 1 — z [57] since at low concentrations there are very few separated 
clusters of spins so that Poo(z) is mostly reduced from unity due to individ- 
ual removed sites. Equation 7, shown as a dotted line in Fig. 8, incorrectly 
predicts p s (z « 30%) = 0, and hence a quantum critical point well below the 
percolation threshold [37]. To our surprise, we find that substituting 1 + z for 
1/P<x>{ z ) quantitatively describes p s (z) / p s (0) even at z = 35%. This modified 
form of Eq. 7 is shown as a solid line in Fig. 8. This substitution matches the 
original expression at low concentrations and prevents the second term in Eq. 
7 from going to zero below the percolation threshold. We note that Eq. 7 is 
a one-loop renormalization-group result, and (unknown) higher-order terms 
might perhaps improve agreement with our observations. Recent numerical 
finite-size scaling results for the spin stiffness give qualitatively similar results 
[57], but deviate from the unmodified form of Eq. 7 more quickly at lower con- 
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centrations. The expression for the spin-wave velocity corresponding to Eq. 7 
is c(z)/c(0) = A(z)(l + z/2) [37], and is shown by the dashed line in Fig. 8. 
This expression does not match our QMC fit results, possibly for a number 
of reasons. Even at z = 0, a fit results in a 20% error in c [13]. Moreover, in 
the presence of random dilution, spin waves become strongly damped and the 
spin-wave velocity is not expected to be a well-defined quantity. 

0. 6 Discussion 

Recent theoretical work for the randomly diluted S — 1/2 NN square-lattice 
Heisenberg antiferromagnet has produced exact results at low concentration 
[33] and led to predictions at intermediate concentrations [37,35]. Impurities 
may localize spin excitations and lead to a breakdown of the classical hydrody- 
namic description of excitations in terms of spin waves above a characteristic 
length scale [35]. Static properties such as the staggered magnetization and 
Tjy, as well as the spin stiffness p s , are expected to remain well defined. On 
the other hand, dynamic observables, such as the spin-wave velocity c, become 
ill-defined in this picture. As long as the ground state remains ordered, the 
low-temperature correlation length should still scale as £ ~ e 27rps / T . This is 
consistent with our observation for £ (z, T) over a wide range of concentrations 
and temperatures, which we find to be well described by Eq. 4, and especially 
the heuristic crossover form Eq. 5. 

Although the percolation transition for the 2D S=l/2 Heisenberg antiferro- 
magnet may be classical [57,58,59], recent QMC simulations of diluted bilayers 
indicate that the percolation threshold for the single layer system is very close 
to a new quantum multi-critical point in an extended parameter space [58,59]. 
Indeed, the dynamic critical exponent obtained in these studies is consistent 
with the thermal correlation length exponent vt ~ 0.7 that we have found 
here. Therefore, in the temperature regime that we have studied, the prop- 
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erties of La 2 Cui_ z (Zn,Mg) z 4 near the percolation threshold appear to be 
controlled by this new critical point. Future inelastic neutron scattering mea- 
surements of the full dynamic structure factor should give further insight into 
this complex quantum many-body system. 

O.P.V. and M.G. thank A. Aharony, A. H. Castro Neto, A. L. Chernyshev, A. 
W. Sandvik, and E. F. Shender for helpful discussions. The work at Stanford 
was supported by the U.S. Department of Energy under contract nos. DE- 
FG03-99ER45773 and DE-AC03-76SF00515, by NSF CAREER Award no. 
DMR9985067, and by the A.P. Sloan Foundation. 

References 

[1] K. B. Lyons, P. A. Fleury, J. P. Remeika, A. S. Cooper, T. J. Negran, Phys. 
Rev. B 37 (1988) 2353. 

[2] Y. Endoh, K. Yamada, R. J. Birgeneau, D. R. Gabbe, H. P. Jennsen, 
M. Kastner, Phys. Rev. B 37 (1988) 7443. 

[3] G. Aeppli, S. M. Hayden, H. A. Mook, Z. Fisk, S.-W. Cheong, D. Rytz, J. P. 
Remeika, G. P. Espinosa, A. S. Cooper, Phys. Rev. Lett. 62 (1989) 2052. 

[4] Y. Tokura, S. Koshihara, T. Arima, H. Takagi, S. Ishibashi, T. Ido, S. Uchida, 
Phys. Rev. B 41 (1990) 11657. 

[5] B. Keimer, N. Belk, R. J. Birgeneau, A. Cassanho, C. Y. Chen, M. Greven, 
M. A. Kastner, A. Aharony, Y. Endoh, R. W. Erwin, G. Shirane, Phys. Rev. B 
46 (1992) 14034. 

[6] M. Greven, R. J. Birgeneau, Y. Endoh, M. A. Kastner, M. Matsuda, G. Shirane, 
Z. Phys. B 96 (1995) 465. 

[7] R. J. Birgeneau, M. Greven, M. A. Kastner, Y. S. Lee, B. O. Wells, Phys. Rev. 
B 59 (1999) 13788. 



16 



[8] Y. J. Kim, R. J. Birgeneau, F. C. Chou, M. Greven, M. A. Kastner, Y. S. Lee, 
B. 0. Wells, A. Aharony, O. Entin-Wohlman, I. Y. Korenblit, A. B. Harris, 
Phys. Rev. B 64 (2001) 024435. 

[9] Y. J. Kim, R. J. Birgeneau, F. C. Chou, R. W. Erwin, M. A. Kastner, Phys. 
Rev. Lett. 86 (2001) 3144. 

[10] R. Coldea, S. M. Hayden, G. Aeppli, T. G. Perring, C. D. Frost, T. E. Mason, 
S. W. Cheong, Z. Fisk, Phys. Rev. Lett. 86 (2001) 5377. 

[11] M. Makivic, H. Q. Ding, Phys. Rev. B 43 (1991) 3562. 

[12] J.-K. Kim, M. Troyer, Phys. Rev. Lett. 80 (1998) 2705. 

[13] B. B. Beard, R. J. Birgeneau, M. Greven, U.-J. Wiese, Phys. Rev. Lett. 80 
(1998) 1742. 

[14] S. Chakravarty, B. I. Halperin, D. R. Nelson, Phys. Rev. B 39 (1989) 2344. 

[15] P. Hasenfratz, F. Niedermayer, Phys. Lett. B 268 (1991) 231. 

[16] S. Tyc, B. I. Halperin, S. Chakravarty, Phys. Rev. Lett 62 (1989) 835. 

[17] R. R. P. Singh, Phys. Rev. B 39 (1989) 9760. 

[18] J. Igarashi, Phys. Rev. B 46 (1992) 10763. 

[19] A. V. Chubukov, S. Sachdev, J. Ye, Phys. Rev. B 49 (1994) 11919. 

[20] C. J. Hamer, Z. Weihong, J. Oitmaa, Phys. Rev. B 50 (1994) 1994. 

[21] N. Elstner, A. Sokol, R. R. P. Singh, M. Greven, R. J. Birgeneau, Phys. Rev. 
Lett. 75 (1995) 938. 

[22] R. R. P. Singh, M. P. Gelfand, Phys. Rev. B 52 (1995) 15695. 

[23] A. Cuccoli, V. Tognetti, R. Vaia, P. Verrucchi, Phys. Rev. B 56 (1997) 14456. 

[24] P. Hasenfratz, Eur. Phys. J. B 13 (2000) 11. 

[25] L. Siurakshina, D. Ihle, R. Hayn, Phys. Rev. B 64 (2001) 104406. 

[26] P. V. Shevchenko, A. W. Sandvik, O. P. Sushkov, Phys. Rev. B 61 (2000) 3475. 

17 



[27] J. F. Annett, R. M. Martin, A. K. McMahan, S. Satpathy, Phys. Rev. B 40 
(1989) 2620. 

[28] A. P. Young, H. Rieger, Phys. Rev. B 53 (1996) 8486. 
[29] S. Sachdev, A. P. Young, Phys. Rev. Lett. 78 (1997) 2220. 
[30] T. Senthil, S. Sachdev, Phys. Rev. Lett. 77 (1996) 5292. 
[31] H.-O. Heuer, Phys. Rev. B 45 (1992) 5691. 

[32] A. W. Sandvik, E. Dagotto, D. J. Scalapino, Phys. Rev. B 56 (1997) 11701. 

[33] S. Sachdev, C. Buragohain, M. Vojta, Science 286 (1999) 2479. 

[34] M. Vojta, C. Buragohain, S. Sachdev, Phys. Rev. B 61 (2000) 15152. 

[35] A. L. Chernyshev, Y. C. Chen, A. H. Castro Neto, Phys. Rev. B 65 (2002) 
104407. 

[36] C. Yasuda, A. Oguchi, J. Phys. Soc. Japan 68 (1999) 2773. 

[37] Y.-C. Chen, A. H. Castro Neto, Phys. Rev. B 61 (2000) R3772. 

[38] A. W. Sandvik, M. Vekic, Phys. Rev. Lett. 74 (1995) 1226. 

[39] R. J. Birgeneau, R. A. Cowley, G. Shirane, H. Yoshizawa, J. Stat. Phys. 34 
(1984) 817, and references therein. 

[40] A. Chakraborty, A. J. Epstein, M. Jarrel, E. M. McCarron, Phys. Rev. B 40 
(1989) 5296. 

[41] S.-W. Cheong, A. S. Cooper, L. W. J. Rupp, B. Batlogg, Phys. Rev. B 44 (1991) 
9739. 

[42] S. T. Ting, P. Pernambuco-Wise, J. E. Crow, E. Manousakis, Phys. Rev. B 46 
(1992) 11772. 

[43] R. L. Lichti, C. Boekema, J. C. Lam, D. W. Cooke, S. F. J. Cox, S. T. Ting, 
J. E. Crow, Physica C 180 (1991) 358. 

18 



[44] G. Cao, J. W. O'Reilly, J. E. Crow, L. R. Testardi, J. Appl. Phys. 75 (1994) 
6595. 

[45] M. Corti, A. Rigamonti, F. Tabak, P. Carretta, F. Licci, L. Raffo, Phys. Rev. 
B 52 (1995) 4226. 

[46] K. Uchinokura, T. Ino, I. Terasaki, I. Tsukada, Physica B 205 (1995) 234. 

[47] M. Hiicker, V. Kataev, J. Pommer, J. Harras, A. Hosni, C. Pflitsch, R. Gross, 
B. Biichner, Phys. Rev. B 59 (1999) R725. 

[48] S. J. Clarke, A. Harrison, J. Mag. Mag. Materials 140 (1995) 1627. 

[49] D. J. Breed, K. Gilijamse, J. W. E. Sterkenburg, A. R. Miedema, J. Appl. Phys. 
41 (1970) 1267. 

[50] D. Stauffer, A. Aharony, Introduction to Percolation Theory, Revised 2 nd Ed., 
Taylor & Francis, Bristol, PA, 1994. 

[51] M. E. J. Newman, R. M. Ziff, Phys. Rev. Lett. 85 (2000) 4104. 

[52] S. Miyashita, J. Behre, S. Yamamoto, in: Computational Approaches in 
Condensed Matter Physics, Springer, Berlin, 1992, p. 97. 

[53] K. Kato, S. Todo, K. Harada, N. Kawashima, S. Miyashita, H. Takayama, Phys. 
Rev. Lett. 84 (2000) 4204. 

[54] S. S. Todo, H. Takayama, N. Kawashima, Phys. Rev. Lett 86 (2001) 3210. 

[55] C. Yasuda, S. Todo, M. Matsumoto, H. Takayama, Phys. Rev. B 64 (2001) 
092405. 

[56] O. P. Vajk, P. K. Mang, M. Greven, P. M. Gehring, J. W. Lynn, Science 295 
(2002) 1691. 

[57] A. W. Sandvik, Phys. Rev. B 66 (2002) 024418. 

[58] A. W. Sandvik, Phys. Rev. Lett. 89 (2002) 177201. 

[59] O. P. Vajk, M. Greven, Phys. Rev. Lett. 89 (2002) 177202. 

19 



[60] A. Aharony, R. J. Birgeneau, A. Coniglio, M. A. Kastner, H. E. Stanley, Phys. 
Rev. Lett. 60 (1988) 1330. 

[61] T. Thio, T. R. Thurston, N. W. Preyer, P. J. Picone, M. A. Kastner, H. P. 
Jennsen, D. R. Gabbe, C. Y. Chen, R. J. Birgeneau, A. Aharony, Phys. Rev. B 
38 (1988) 905. 

[62] A. P. Murani, A. Heidemann, Phys. Rev. Lett. 41 (1978) 1402. 

[63] J. D. Axe, A. H. Moudden, D. Hohlwein, D. E. Cox, K. M. Mohanty, A. R. 
Moodenbaugh, Y. Xu, Phys. Rev. Lett. 62 (1989) 2751. 

[64] M. K. Crawford, R. L. Harlow, E. M. McCarron, W. E. Farneth, J. D. Axe, 
H. Chou, Q. Huang, Phys. Rev. B 44 (1991) 7749. 

[65] H. G. Evertz, G. Lana, M. Marcu, Phys. Rev. Lett. 70 (1993) 875. 

[66] U.-J. Wiese, H.-P. Ying, Z. Phys. B 93 (1994) 147. 

[67] M. Greven, R. J. Birgeneau, Phys. Rev. Lett. 81 (1998) 1945. 

[68] E. Manousakis, Phys. Rev. B 45 (1992) 7570. 

[69] A. H. Castro Neto, D. Hone, Phys. Rev. Lett. 76 (1996) 2165. 

[70] R. J. Birgeneau, R. A. Cowley, G. Shirane, H. J. Guggenheim, Phys. Rev. Lett. 
37 (1976) 940. 

[71] A. B. Harris, S. Kirkpatrick, Phys. Rev. B 16 (1977) 542. 



20 



Figure Captions 



Fig. 1: Schematic of finite-sized sections of the infinite square lattice with random site dilution. Sites 



Fig. 2: Measurements of the Neel temperature Tm(z) for various samples, a) Magnetic susceptibility r 



Fig. 3: Order parameter measurements performed on the 35% sample on the cold neutron instrument 



Fig. 4: Neutron diffraction measurements of the first-order structural phase transition from the orthoi 



Fig. 5: Magnetic and structural phase diagram of La9Cui_.,fZn.Mg^0^. a) Neel temperature vs. dilut 



Fig. 6: Two-axis, energy-integrating measurements of the static structure factor taken with 30.5 meV 



Fig. 7: a) Spin-spin correlation lengths, in units of the lattice constant, versus inverse temperature, in 



Fig. 8: Effective spin-wave velocity c(z) and spin stiffness 2ttoJz) as a function of non-magnetic conc( 
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